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1 Introduction 



Polygonal Markov fields, originally introduced by Arak & Surgailis [1, 2] and further stud- 
ied by Arak, Clifford & Surgailis [3, 4] are random ensembles of non-intersecting polygonal 
contours in the plane interacting by hard-core exclusions and exhibiting two-dimensional 
germ-Markov property [2], with a variety of additional possible terms entering the Hamil- 
tonian, including length and area elements (ibidem). The polygonal fields with V-shaped 
nodes (no nodes of order higher than two) as considered in this paper, share a number 
of essential features with the two-dimensional Ising model, prominent examples including 
the presence of an Ising-like phase transition [15, 16] as well as low temperature phase 
separation and Wulff droplet creation [17]. For these reasons, the polygonal Markov fields 
are sometimes regarded as continuum counterparts of the Ising model in the plane (as 
well as of the Potts model if higher order nodes are admitted). Remarkably, in many as- 
pects the polygonal fields are exactly tractable, especially in the so-called consistent regime 
falling into the supercritical temperature region. In particular, at the consistency point 
one knows the exact value of the partition function, first and second order characteristics 
of the field [2, 3, 4]. Further, a lot is known about the higher order correlations as well, 
including certain exact formulae [18] and necessary and sufficient conditions for factorisa- 
tion of edge correlations (ibidem) as well as an exponential mixing statement (asymptotic 
factorisation) for edge correlations (ibidem). A striking feature of polygonal Markov fields 
is that they admit a number of particularly convenient algorithmic constructions - graph- 
ical representations [2, 3, 4, 12, 16, 17, 18] which are in fact the main tool for establishing 
of the afore-mentioned results. The geometric ingredient in these considerations is so pre- 
dominant that often no supplementary calculations are needed, which stands in a strong 
contrast to the classical Ising model. 

The rich class of graphical constructions developed for polygonal fields have also found 
their applications in Bayesian image processing where they are used to generate image 
segmentations, see [6, 10, 11, 12, 19]. Experimenting with various black-white and grayscale 
images we already obtained promising results, further algorithms are a subject of our 
ongoing research in progress. 

The purpose of the present paper is to complement the existing exact results for con- 
sistent polygonal Markov fields by establishing in Theorem 2 an explicit stochastic repre- 
sentation for their higher order edge correlations in terms of expectations of the so-called 
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crop functionals of polygonal webs. The polygonal web arises as the union of a collection 
of trajectories of continuous time critical branching polygonal random walks in the plane, 
interacting by a bridge- creating mechanism attempting to clasp random trees generated 
by branching walks into a web by establishing linear bridges between trees. As mentioned 
above, the polygonal fields admit a graphical construction whose full details in Section 3 
below. The dynamics of this construction can be interpreted as a (quite untypical) inter- 
acting particle system evolving in time and, speaking in this vein, the polygonal web arises 
in a graphical construction, given in Section 4 below, which can be regarded as dual to 
that for polygonal fields. The nature of this duality consists, roughly speaking, in the fact 
that the dynamic representation of polygonal webs exhibits features strongly reminiscent of 
those of the polygonal field construction under inverted time flow direction. To some extent 
this can be perceived as an analogy to the classical duality for interacting particle systems, 
see e.g. Section III. 4 in [13], although this is not a close analogy, expecially that here we 
deal with entire histories of the considered interacting particle systems (whose trajectories 
trace the polygonal fields and webs) rather than just with their instantaneous configura- 
tions. The proof of our duality representation for edge correlations goes by constructing 
a martingale interpolating between edge correlation functions and crop functionals of the 
corresponding polygonal webs. 

The rest of this paper is organised as follows. In Section 2 below we present the 
concept of a non-homogeneous polygonal field as introduced in [18]. Next, in Section 3 we 
discuss the generalised graphical construction of such fields as developed ibidem. In the 
further Section 4 we develop a dual graphical construction and define the polygonal web. 
In the next Section 5 we define edge correlation functions of polygonal fields and state 
our main representation Theorem 2 as briefly discussed above. The proof of this theorem 
is given in Section 6, where a number of auxiliary constructions are also developed and 
many auxiliary results established. Finally, to the last Section 7 we delagate the proof of 
a technical existence result for edge correlations. 

2 Non-homogeneous polygonal Markov fields in the 
plane and their consistent regime 

For an open bounded convex set D define the family T D of admissible polygonal configura- 
tions on D by taking all the finite planar graphs 7 in D U dD, with straight-line segments 
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as edges, such that 

• the edges of 7 do not intersect, 

• all the interior vertices of 7 (lying in D) are of degree 2, 

• all the boundary vertices of 7 (lying in dD) are of degree 1, 

• no two edges of 7 are colinear. 

In other words, 7 consists of a finite number of disjoint polygons, possibly nested and 
chopped off by the boundary. Further, for a finite collection (/) = (k)f =1 of straight lines 
intersecting D we write T D {1) to denote the family of admissible configurations 7 with the 
additional properties that 7 C lj™ =1 U and 7 fl U is a single interval of a strictly positive 
length for each = 1, ...,n, possibly with some isolated points added. 

For a Borel subset of A C R 2 by [[A]] we shall denote the family of all straight lines 
hitting A so that in particular [[R 2 ]] stands for the collection of all straight lines in R 2 . We 
shall also write [[A]] for the family of all linear segments in R 2 hitting A. Further, let p 
be the standard isometry-invariant Haar-Lebesgue measure on the space [[R 2 ]] of straight 
lines in R 2 . Recall that one possible construction of p goes by identifying a straight line / 
with the pair (0, p) e [0, n) x R, where (p sin(0), p cos(0)) is the vector orthogonal to /, and 
joining it to the origin, and then by endowing the parameter space [0, ir) x R with the usual 
Lebesgue measure. Note that the above parametrisation of [[R 2 ]] with [0, tt) x R endows 
[[R 2 ]] with a natural metric, topology and Borel cr-field which will be used in this paper. 
Next, consider a non-negative Borel measure Ai on [[R 2 ]] admitting a locally bounded 
density m(-) with respect to p. Below, the measure M. will be interpreted as the activity 
measure on [[R 2 ]]. Let A M be the Poisson line process on [[R 2 ]] with intensity measure 
M. and write A^ 1 for its restriction to the domain D. Further, define the Hamiltonian 
L M : Yd — > R + given by 

L*(7):= E ^M). ^ er ^- (!) 

eG Edges (7) 

We note that the energy function L M should be regarded as an anisotropic environment- 
specific version of the length functional. Indeed, interpreting the activity Ai(dl) of a line 
/ hitting an edge e G Edges(7) as the likelihood of a new edge being created along / 
intersecting and hence fracturing the edge e in 7, we observe that, roughly speaking, the 
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value of -M([[e]]) determines how likely the edge e is to be fractured by another edge present 
in the environment. In other words, L M (j) determines how difficult it is to maintain the 
whole graph 7 G without fractures in the environment M. - note that due to the 
anisotropy of the environment there may be graphs of a higher (lower) total edge length 
than 7 and yet of lower (higher) energy and thus easier (more difficult) to maintain and 
to keep unfractured due to the lack (presence) of high activity lines likely to fracture their 
edges. 

Following [18], with the above notation, for j3 G M further referred to as the inverse 
temperature (from mathematical viewpoint also the unphysical negative values of inverse 
temperatures are admissible), we define the polygonal field A^' 13 in D with activity mea- 
sure M. to be the Gibbsian modification of the process induced on by A^ 1 , with the 
Hamiltonian L M at inverse temperature f3, that is to say 

X M-,^M ._ E E 7e r c (A-)nc-p(-^( 7 )) 



P(A D ' e G) :- exp ( _^ (7)) ( 2 ) 

for all sets G C Yd Borel measurable with respect to, say, the usual Hausdorff distance 
topology. Note that this definition can be rewritten as 

P(< ;/3 G d 7 ) oc exp(-^( 7 )) J] M{dl[e\), 7 G T D , (3) 

e£Edgcs(7) 

where l[e] is the straight line extending e. In other words, the probability of having A^' 13 G 
d'-f is proportional to the Boltzmann factor exp(— f3L M {^)) times the product of edge 
activities Ai(dl[e]), e G Edges( 7 ). Observe also that this construction should be regarded 
as a specific version of the general polygonal model given in (2.11) of [2]. The finiteness of 
the partition function 

Z£*:=E Yl exp(-/3L^( 7 ))<oo (4) 
7 er D (A^) 

for all (3 G M is not difficult to verify and has been established in [18], see (32) there. 

In this paper we shall focus on polygonal fields in their consistent regime corresponding 
to /3 — 1. As shown in Section 3 in [18], this particular choice of temperature parameter 
places us in the context of a non-homogeneous version of Arak-Surgailis [2] construction 
for the consistent polygonal fields, see Section 4 there. This ensures striking properties of 
the field. First of these, the consistency, states that for each open bounded and convex 
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DCD'cR 2 the field A% := A^'' 1 coincides in law with A%, : = A^,' 1 n D, thus allowing 
for a direct construction of the infinite volume process (thermodynamic limit) A M := A M '^ 
on the whole M 2 such that Aj? = A M (ID. The infinite volume process A M takes its values 
in the space Y := T K 2 of whole-plane admissible configurations, with obvious meaning of 
this notation. Further, the explicit formula for the partition function is known for 
convex D, see Theorem 4.1 in [2] for the homogeneous case and Theorem 1 in [18] for the 
general non- homogeneous set-up. We state this formula in (6). Moreover, one-dimensional 
linear sections of the field are fully characterised in distribution, see ibidem and Theorem 

2 in [18]. Finally, the polygonal fields Ap enjoy the two-dimensional Markov property 
stating that the conditional behaviour of the field inside a smooth closed curve depends 
on the outside configuration only through arbitrarily small neighbourhoods of the curve 
or, equivalently, through the trace of the external configuration on the curve, consisting of 
intersection points and directions. This property is a direct consequence of the Gibbsian 
definition (2,3) of the field and, unlike the previous properties, it holds for all inverse 
temperatures f3 e R rather than just for (3 = 1. We will not discuss this property in the 
present paper and we refer the reader to the original work of Arak and Surgailis [2] for 
further details. 

3 Generalised dynamic representation for consistent 
fields 

The present section is meant to recall the generalised dynamic representation for consistent 
polygonal fields as developed in Section 4 of [18], which will serve as a crucial tool in our 
further considerations. The name generalised representation comes from the fact that it 
generalises the original construction of such fields introduced by Arak and Surgailis in [2] . 
In the sequel we will often omit the qualifier generalised for the sake of terminological 
brevity. To describe the representation, fix the convex field domain D and let (D t ) te [ 0) i] be 
a time-indexed increasing family of compact convex subsets of D, eventually covering the 
entire D and interpreted as a growing window gradually revealing increasing portions of 
the polygonal field under construction in the course of the time flow. In other words, under 
this interpretation, the portion of a polygonal field in a bounded open convex domain D 
uncovered by time t is precisely its intersection with D t . To put it in formal terms, consider 
(A)te[o,i] satisfying 
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(Dl) (D t ) t6 [ 0; i] is a strictly increasing family of compact convex subsets of D = D U dD. 
(D2) D is a single point x in D = D U <9L>. 
(D3) Di coincides with D. 

(D4) D t is continuous in the usual Hausdorff metric on compacts. 

Note that the extra fifth condition imposed on D t in Section 4 of [18] is automatically 
satisfied here due to the absolute continuity M. <C fi and thus is not mentioned here. 
Clearly, under these conditions, for .M-almost each I G [[D]] the intersection lC\D n consists 
of precisely one point A(Z), where 77 = inf{t G [0, 1], D t C\ I ^ 0}. The point A(/) will be 
referred to as the anchor point for I, this induces the anchor mapping A : [[D]] — >• D defined 
.M-almost everywhere. Consider now the following dynamics in time t G [0,1], with all 
updates given by the rules below performed independently of each other. 

(GE:Initialise) Begin with empty field at time 0, 

(GE:Unfold) Between critical moments listed below, during the time interval [t,t + dt] 
the unfolding field edges in D t reaching dD t extend straight to D t+dt \ D t , 

(GE:BoundaryHit) When a field edge hits the boundary dD, it stops growing in this 
direction (note that .M-almost everywhere the intersection of a line with dD consists 
of at most two points), 

(GE:Collision) When two unfolding field edges intersect in D t+dt \ D t , they are not 
extended any further beyond the intersection point (stop growing in the direction 
marked by the interesction point), 

(GE:DirectionalUpdate) A field edge extending along / G [[D t ]] updates its direction 
during [t,t + dt] and starts unfolding along /' G [[Z'*'*" 1 "^]], extending away from 
the anchor point A(/'), with probability M(dl'), where ft' t+dt ^ := I n (D t+dt \ D t ). 
Directional updates of this type are all performed independently, 

(GE:LineBirth) Whenever the anchor point A(Z) of a line / falls into D t +dt \ D t , the 
line / is born at the time t at its anchor point with probability Ai(dl), whereupon it 
begins extending in both directions with the growth of D t (recall that / is .M-almost 
always tangential to dD t here), 
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(GE:VertexBirth) For each intersection point of lines h and / 2 falling into D t+dt \ D t , 
the pair of field lines li and l 2 is born at l\ n l 2 with probability Ai(dli)Ai(dl 2 ), 
whereupon both lines begin unfolding in the directions away from their respective 
anchor points A(li) and A(Z 2 ). 

Observe that the evolution rule (GE:VertexBirth) means that pairs of lines are born at 
birth sites distributed according to a Poisson point process in D with intensity measure 
given by the intersection measure ((M)) of M. defined as follows 

((M))(A) := X -M x M({(h,k), hnl 2 c A}). (5) 

The importance of the intersection measure lies in the fact that 

Z^ = eM((M))(D)) (6) 

as shown in Theorem 1 in [18] . The following result stating that the polygonal field resulting 
from the above construction actually coincides with Aj) has been established in [18], see 
Theorem 3 there. 

Theorem 1 The random contour ensemble resulting from the above construction (GE) 
coincides in law with A^. 

4 Polygonal web 

Having defined the non-homogeneous polygonal fields and presented their graphical con- 
struction, we pass now to another object central to this paper, which we name the polygonal 
web. Whereas the details of the connection between the critical polygonal web and the 
corresponding consistent polygonal field are to be established in the subsequent Section 5, 
here we emphasise that, in a sense, the polygonal web constitutes the dual object to the 
polygonal field sharing the same activity measure, and this duality is going to be reflected 
in the construction of the polygonal web. Roughly speaking, the polygonal web arises 
as the union of critical branching polygonal random walks, interacting by an additional 
bridge- creating mechanism, clasping the branched polygonal trees into a web. 

4.1 Constructing the polygonal web 

Consider an open bounded and convex domain D and let (D t ) te \ 0t i\ be a growing family 
of compact subsets of D satisfying the usual conditions (Dl-4) as in Section 3. Further, 
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assume a collection (l,x) = (k,Xi)f =l is given where Xi are points in D whereas U are 
straight lines with x^ £ k, i — 1, . . . , k. Moreover, to avoid uninteresting pathologies we 
require that no three different lines U, i = 1, . . . , k, intersect at one point. Each such pair 
(li, Xi) will be called an edge marker because k can be interpreted as a directional indicator 
for a linear segment /edge passing through x iy see Section 5 below where this interpretation 
is further developed and exploited. The entire collection (l,x) will be referred to as the 
edge marker configuration. The polygonal web W[(Z,x)] := W^[(/,x)] generated by (l,x) 
in D, with activity measure Ai, is the union of polygonal trees in D arising as the final 
state Wi of the following graphical construction process w s evolving for s £ [0, 1], where all 
random updates listed are performed independently. 

(W: Start) At the time s = we let w consist of zero- length edges (edge germs) at 
Xi, 1 = 1,... ,k, directed along the respective /j's. 

(W:GrowInwards) Between the critical moments listed below, during the time interval 
[s, s + ds] all edges of w s reaching the boundary dD\^ s extend straight to Di- S \ 
Dis-ds along their respective directional lines. The edges (edge germs) not yet 
touched by the boundary of the shrinking domain Di_ s (and hence contained in the 
interior of Di- S ) remain intact and do not evolve until eventually hit by the boundary 
at some later time (unless they get frozen prior to that, see below, in which case they 
never start evolving). Below, we call edges reaching the current boundary dDi_ s 
active and we say that edges (edge germs) not yet hit by the boundary are inactive. 
Note that the intersection point of a web edge with the current boundary dD\^ s can 
be interpreted as its instantaneous growth-tip and so will it be called in the sequel. 
Observe that inactive edges do not have growth-tips. 

(W: Branch AndTurn) During the time interval [s, s + ds] an active web edge reaching 
the boundary dDi_ s and extending along I £ [[£)i_ s ]] yields a new offspring edge 
starting at lt s > s + ds ^ ■= I n D x _ s \ and directed along V £ [[ll s ^+ ds }]] w j t h 

probability M(dl'). Both the original and offspring edges go on evolving according 
to the usual rules. The branching updates are performed independently for all active 
web edges. 

(W:ForcedBranch AndTurn) If during the time interval [s, s + ds] an active web edge 
extending along some I £ [[Z)!_ s ]] and reaching the boundary dDi_ s intersects the 
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directional line I' of some other non-frozen web edge in w s fl Di- S , be it active or 
inactive, and there is no web edge along /' reaching I fl dDi_ s and created in a prior 
stage of the evolution, then a new offspring edge is created at IddD^s directed along 
I 1 and both the original and offspring edges go on evolving according to the usual 
rules. Note that we only consider directional lines I' of edges hitting the domain Di_ s 
- the edges present in D \ Di_ s but terminated before the time s (see below for edge 
termination events) are not taken into account. Thus, the forced branching occurs if 
the growth-tip of an edge hits the current directional line of another (non-frozen) edge 
currently present in the system. Whenever a forced branching occurs, the resulting 
offspring edge is called a forced edge whereas the edge whose directional line /' gave 
rise to the branching is referred to as the corresponding forcing edge. Observe that a 
pair of forced and forcing edges will meet and coalesce into a single edge at further 
stages of the construction unless one of the edges terminates prior to that. 

Note that if the forcing edge is active, it resides at the boundary of the domain D s 
and yields a single forced edge on the opposite side of the domain. On the other 
hand, an inactive forcing edge (edge germ) located in the interior of D s may give rise 
to (at most) two forced edges, one on each side of the domain. 

(W: Terminate) During the time interval [s, s + ds] an active web edge reaching the 
boundary dD\^ s and extending along I G [[Z)i_ s ]] terminates (stops evolving) with 
probability M([[l [s ' s+ds] ]]) where li s ' s + ds ^ ■= I n D^ s \ L>i_ s _ ds . 

(W:StopIfSeparated) Whenever at some time moment s G [0, 1] a web edge (or edge 
germ) e in w s nDi_ s , be it active or inactive, has the property that /[e]nconv([w s \e]fl 
dDx_ s ) = (that is to say the directional line l[e] of the web edge e does not hit the 
convex hull generated by the current growth-tips of the remaining non-terminated 
active edges and non-frozen germination points of the remaining inactive edges, in 
which case we say that e separates from w s at the time s), then e terminates and 
stops evolving at this point. Note that in case of an inactive edge germ e to terminate 
means to remain frozen in inactive state and never to activate even when hit by the 
boundary at the later stages of the evolution. 

A careful reader might ask at this moment why in the above construction we do not consider 
the case when at some time s a web edge becomes tangential to the domain Di_ s . The an- 
swer is that, with probability one, such cases do not occur in the course of the construction 
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because an edge to become tangential to the boundary of the domain separates from the 
web prior to that and thus gets terminated by an application of (W:StopIfSeparated) 
rule. 

The construction of the polygonal web as presented above admits a natural intuitive 
description - the edge germs initiating the process emit critical branching polygonal ran- 
dom walks directed by the activity measure M. and unfolding inward the domain D which 
can be regarded as dual to the dynamic representation in Section 3 where the growth was 
directed outwards. The branching is critical because the binary branching and termina- 
tion intensities coincide in (W:BranchAndTurn) and (W:Terminate). The role of the 
additional (W:ForcedBranchAndTurn) rule is to ensure the possibility of bridging the 
gaps between two separated co-linear parts of the same segment present on two opposite 
sides of the window Di_ s , thus clasping the polygonal trees into a web. Finally, as may be 
seen in the sequel, the (W:StopIfSeparated) rule reflects the structural knowledge on 
independence of edge covering events in polygonal fields as established in [18], and as such 
it is not indispensible in its full form for the theory developed below to be valid and may 
be replaced by various weaker variants, see Remark 1. 

It is useful to note that, regarded as a polygonal graph, the polygonal web contains 
T-shaped nodes (branching points), I-shaped nodes (edge terminal points) and X-shaped 
nodes (edge intersection points) but no V-shaped nodes. 

4.2 Crop functional 

To establish a direct link between the polygonal web and polygonal fields we define now 
the crop functional of the polygonal web, further denoted as crop(W[(I, x)]). To this end, 
we identify the polygonal web W[(Z, x)] = yV^[(J, x)\ with the collection of branches con- 
necting the initial edge germ locations Xi, % — 1,... ,k to the edge termination points 
Vi, ■ ■ ■ ,Um resulting from (W:Terminate,StopIfSeparated) rules or arising as meeting 
points where forcer-forced pairs of co-linear web edges unfolding in opposite directions 
merge into linear segments. To avoid nuisance technicalities below we formally interpret 
each meeting point of a forcer-forced pair as two distinct points, one terminating the forcing 
branch, the second terminating the forced branch. Clearly, m > k since each initial edge 
germ x^ emits at least one branch and each branch eventually terminates. Keeping (l,x) 
fixed, we shall index the branches constituting W[(/,a;)] by their terminal points, writing 
branch [yj] for the branch terminating at i/j - the inambiguity of this indexation is ensured 
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by the above convention on meeting points of forcer-forced pairs. Along each branch we 
have a natural chronological ordering from the root x iy i G {1, ... ,k} to the endpoint 
Uj, j G { 1 , • • • , m}. For a collection of branches branch[y], y <E y, y C . . . ,y m }, we 
consider the induced polygonal crop graph CropGraphf^] obtained as follows. 

(Crop:Grow) Follow the growth of all branches in {branch[y], y G 3^} starting from 
their roots and unfolding towards their respective endpoints during the time interval 
[0, 1] as in the course of the polygonal web dynamics (W). 

(Crop:StopOnCollision) Whenever in the course of their growth two branches meet, 
they both stop growing at this point. 

Thus, the crop graph is a subgraph of W[(l, x)] containing T-shaped, I-shaped and V- 
shaped nodes but not X-shaped nodes. This is because the crop graph arises by (recur- 
sively) cutting off the parts of branches past their intersections with other branches present 
in the inducing collection. There are two ways in which two branches can meet in the above 
construction - they can either intersect coming from two non-colinear directions and giving 
rise to a V-shaped node in CropGraphf^], or meet coming from opposite co-linear direc- 
tions yielding a linear segment rather than a graph node. Note that two distinct branches 
sharing a common sub-branch and thus coinciding during initial growth phase are not 
considered to meet or intersect! On the other hand, if several distinct branches coinciding 
during an initial growth phase intersect another branch(es) during this phase, the growth- 
interrupting (Crop:StopOnCollision) rule applies to all these branches simultaneously. 

It is easily seen that two different collections of terminal points can yield identical crop 
graphs because some branch turning points can be cut off due to collisions in the course 
of the (Crop) dynamics. It is clear though that for each instance of a crop graph arising 
in (Crop) dynamics there exists a unique collection of terminal points with the property 
that no branch turning points occur past the cut-off points. This unique collection is called 
minimal for its crop graph, or just minimal for short if no ambiguity arises. 

To proceed, we say that a collection y C . . . , y m } of branch-determining endpoints 
is complete iff each initial germ location is the root of some branch [y], y G y. Clearly, the 
cardinality of such a complete collection cannot fall below k. Further, we say that the crop 
graph CropGraphf^] of a complete endpoint collection y is normal iff it contains no forced 
edges which fail to eventually meet and merge with their co-linear forcing edges. Note that 
this condition can be violated by either having a forced edge without its forcer present in 



12 



the graph or due to a death or directional update along either a forcing or a forced edge. 
Whereas the death of either edge in an forcer-forced pair does necessarily lead to the lack 
of normality, a directional update does so only if there is no other branch in the collection 
along which the considered forcer-forced pair could extend further past the turning point. 
A graph which is not normal is called abnormal. Thus, roughly speaking, an abnormal 
graph is a graph containing the forced part of an incomplete bridge. Write 

, , J 1, if y C {y 1 , . . . ,y m } is complete, minimal and CropGraphf^] is normal, 
^ ' \ 0, otherwise. ^ 

With this notation, we define the crop functional 

crop(W[(/»]) := {-l) caTd(y) - k t{y)- (8) 

y>Q{yi,--- ,ym} 

The expression (8) has the aesthetic advantage of defining the crop functional in a form 
reminiscent of the classical inverse Mobius transform, with the summation performed over 
all subsets of {y±, . . . ,y m }, see e.g. Section 2.6 in [14]. To exploit this feature define 



1, if CropGraph[y] is normal, 
0, otherwise, 



that is to say i is the indicator of crop graph normality, without the extra completeness and 
minimality requirements. Given a crop graph g = CropGraphf^] for some complete and 
minimal y e C {y 1 , . . . , y m }, we let B[y], y G y g , be the set of all yj, j = 1, . . . , m, such that 
branch^] contains the entire subbranch branch [|/]n^. Then, for any y C [y ly . . . , y m ] such 
that CropGraph[y] = CropGraph^J we have i(y) = t(y e ) and, moreover, y decomposes 
into the disjoint union of non-empty y y = y fl B[y], y E y e . We also have card(y) = 
card(3^) + J2 y ey e l CSiT( ^0^y) ~~ !]■ Consequently, using Newton's binomial formula, 

j2 (-i) ca * d w- k t(y) = 

y, CropGraph[y]=CropGraph[ye] 
(_ 1 )card(y e )-fc t (-y^ "Q ^ ( _ 1)C ard(^)-l = (_ l^TO-fc^) "Q j = 

yey e $^y y CB[ y ] y&y e 

(-i) caid ^ k L(y e ). 

Thus, (8) can be alternatively rewritten as 

crop(W[(Z»])= (-l) caid{y) - k t(y). (9) 

yc{ yi ,..., ym },y complete 
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Further, we put T[xj\ := {yj G {yi, . . . ,y m }, root (branch [y 3 -]) = Xi}, i = 1,... ,k, and 
write, applying Newton's binomial formula, 

k k 

^ ^_ 1 ^card(y)-fe = JJ Y (_iyzrd(yi)-l = J~£ 1 = 1. 

ycfo,.. >Wm }, y complete i=1 fc^cr^] *=i 

Combining this with (9) we conclude that 

crop(W[(/»]) = 1- ^ (-l) cajd ^- fc [^) - 1] (10) 

yc{ yi ,..., ym },y complete 

and hence the crop equals one for polygonal webs whose all complete branch subcollections 
yield normal crop graphs, whereas the deviations of the crop functional from the standard 
value one are due to crop graph abnormalities, which will be further exploited in the sequel. 

A natural alternative way of defining the crop functional involves summation over crop 
subgraphs of W[(l, x)], that is to say over all possible different graphs arising in the (Crop) 
dynamics above, in which case (8) becomes 

crop(W[(I, x)\)= nUmb6r ° f branchingS in e 

e is a normal crop graph in w[(i (11) 

since the number of branchings in g = CropGraph[3^], y minimal, is easily seen to coincide 
with card(^y) — k. 

4.3 Interpretation of the crop functional 

A few words are due at this point to provide an intuitive interpretation of the crop func- 
tional as formally defined in (8). To this end, we begin by mentioning a kid game quite 
popular in the happy time of the author's childhood: given a collection of arrows on a 
sheet of paper draw a family of closed curves passing through these arrows, and in case 
where this can be done in more than one way resolve the ambiguity by trying to make 
the resulting picture resemble some real-life object. In fact, this and related problems are 
not just games and find serious interest in studies on human and computer vision, see 
[5] and the references therein. In mathematical terms and specialising to our polygonal 
set-up, given an edge marker collection (l,x) = (k,Xi)f =l we ask for admissible polygonal 
configurations 7 G T R 2 with the following properties 

• Each edge e of 7 contains some edge marker point x^ e ) such that l[e] = k( e )- 
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• Each edge marker point Xi is contained in some e(i) G Edges(7) such that l[e(i)] = k. 



Following Section 5 of [18] we denote the family of such configurations by T(l,x) and 
write N(l,x) for the cardinality of F(l,x), that is to say the number of solutions to the 
discussed problem. To proceed, consider first a particularly simple deterministic instance 
W°[(Z,x)] := W£)[(Z,x)] of polygonal web generated by (l,x) not depending on the activity 
measure M. - define W°[(Z, x)\ to arise in the course of the above (W) dynamics without 
non-forced turns/branchings and without termination events, which is a usual situation 
for example when x^s are very close to each other and are all contained in a domain D 
with very small -M([[.D]]) where applications of (W:TurnAndBranch, Terminate) are 
very unlikely. Observe that the notation comes from the fact that the above 

construction of this polygonal web coincides with the (W) dynamics under zero activity 
measure. Then, as will be shown in Lemma 2 below, 

crop(W° [(l,x)}) =N(l,x) 

for each bounded convex domain D containing all x^s. 

At the other extremity lies the typical behaviour of the crop functional of polygonal 
webs generated by marker configurations (l,x) where the distances between all X{ s are 
huge and where all Zj's are pairwise different - due to the criticality of our branching 
mechanism the usual situation then is that the web W[(Z, x)] splits into disjoint and distant 
sub- webs originating from respective x^s and the bridge creating attempts between these 
sub- webs fail with overwhelming probability in consequence of (W:Terminate). In these 
circumstances all normal crop graphs in W[(Z,x)] are readily seen to be unions of disjoint 
normal crop graphs in sub-webs stemming from individual Xj's and thus, by (11), the value 
of crop(W[(Z, x)]) factorises into the product of crops of the sub-webs. Using (15) below 
combined with our main Theorem 2 and resorting to Remark 1 allows us to conclude 
easily but not immediately that the expectations of these individual crops are all 1 and 
consequently 

Ecrop(W[(Z,x)]) « 1 

as well. This fact, not studied in detail here, is intimately related to mixing properties 
of polygonal fields (asymptotic factorisation of edge correlations) but this link will not 
be followed any further in this paper because so far we are only able to establish a slow 
polynomial mixing using a rather technical argument based on polygonal webs, whereas 
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alternative methods developed in Section 7 of [18] allowed us to establish exponential mixing 
there at least for rectangular fields. 

In remaining situations the crop functional interpolates between the above extremities. 

5 Polygonal web representation for edge correlations 

The purpose of this section is to formulate our main result stating that arbitrary order edge 
correlation functions of the polygonal field coincide with the expectations of the respective 
polygonal web crops. 

Edge correlations Having introduced the crucial concepts in preceding sections we are 
now in a position to define the principal object of our study in this paper, that is to say 
the edge-correlation functions for polygonal fields, and to formulate our main results. Due 
to the polygonal nature of the considered field the natural object to consider are the edge 
correlations 

a M [dl 1 ,x 1 ; . . . ;dl k ,x k ] : = P (Vf =1 3 eeEdgcs( . 4 .M ) ir h (xi) G e, l[e] G dl t ) , (12) 

where h, . . . ,l k are straight lines and 7T; i is the orthogonal projection on U. In all cases 
below we shall be interested in correlations with Xi G h, in which case a M [dl±, xi, . . . ; 
dlk,Xk] can be interpreted as the probability element that the polygonal field A M passes 
through points Xi in the directions determined by the respective lines l iy i = l,...,k. 
For general Xi, not necessarily lying on / i? the A;-fold correlation a M [dh, x±] . . . ;dlk,Xk] 
is the probability that the polygonal field passes through points ttiX x i) m the directions 
determined by the respective lines i — 1, . . . , k. 

Recall from Section 4 above that collections (T, x) = a; i )f : =1 of lines k and points Xi 
with Xi G U are referred to as edge marker configurations (collections) whereas each pair 
(h,Xi) belonging to such a collection is called an edge marker. No edge marker can occur 
twice in an edge marker collection. An edge along li passing through xi is said to cover 
the marker We say that an edge marker collection (I, x) = (xi,h)f = i is in general 

position if the lines U are pairwise different and Xj G" U for j ^ i, otherwise if U = lj for 
some i ^ j then the collection is called degenerate and declared coupledhy k = lj, 

finally if X{ G lj for some i ^ j with /j ^ lj then the collection is said to be in singular 
position. Thus, an edge marker collection can be simultaneously degenerate and singular. 
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As mentioned above, if two edge marker lines U = lj in a degenerate configuration coincide, 
we say that the edge markers (li,Xi) and (lj,Xj) are coupled, sometimes for brevity we just 
say that X{ and Xj are coupled. While allowing both for singularity and degeneracy of 
marker collections, we strictly exclude the situations where three or more different marker 
lines intersect at one point, in order to avoid unnecessary technical pathologies. 

Note that in the singular case where Xi G lj for some i ^ j but Zj ^ lj it makes 
sense to consider one-sided correlations a M [dlf, x^, dlj, Xj, . . . ] and a M [dl~ , xf, dlj, Xj, . . . ] 
with If and l~ standing for two half-lines into which /j is split by the intersection point 
{x^ — Zj fl lj. The definition of such correlations is analogous to (12) the difference being 
that the field edge containing Xj is required to extend respectively at least in the direction 
of If and l~ , yet it is also allowed although not required to extend in the opposite direction 
as well. Since M. <C /x, it is easily seen that for a configuration in general position we would 
have a [dlf, xf, dlj, Xj] . . . ] = cr M [dl~ ,xf, dlj, Xj] . . .] = a M [dk, xf, dlj, Xj\ . . . ], but this is 
no more the case in the considered singular situation, where the respective correlations are 
non-trivially affected by the event that an edge along lj may extend from Xj down to x^ 
where it may intersect with another edge along /j. 

In the sequel we will use the notation A* [A;], * e {g,s,d}, for the respective sets of 
all collections (/j,Xj)f =1 in general, singular and degenerate positions. Further, we put 
A* := IXi^m * G {g,s,d} and A[k] := \J* e {g,s,d} A *M ^d A := [j^ s4} A[k]. 
Clearly, A[k] can be endowed with the natural product topology from ([[1R 2 ]] x R 2 ) xfc . 

Edge correlation functions For a non-degenerate edge marker configuration (I, x) = 
{hi x i)i=i the (normalised) k-fold edge correlation function (j>(li,xi, ... ; 4, Xfc) is defined by 

<P(h, Xl ; ...;l k ,x k ) = ^(h,x 1 ;...; l k , x k ) := °- : _1 ^ip! , (13) 

More generally, if (I, x) is a degenerate edge marker configuration, its normalised correlation 
is given by 

Mf ^ _ o- M [dli,xu... ;dl k ,x k ] 

YTM(dh) l j 

where Yi* stands for the product over lines U in which each line is present exactly once, 
with repetitions discarded. We also adopt the convention that 0(0) := 1, where stands for 
empty edge marker configuration. The existence of edge correlation functions is guaranteed 
by the following lemma. 
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Lemma 1 The correlation function <j)(l, x) exists for each edge marker configuration (x, I) 
and is continuous on A 9 [k]. In particular, <f>(-) is locally bounded. 

Likewise, we can also consider one-sided versions of the correlation functions, for which 
the existence and local boundedness statements of Lemma 1 readily extend. Note that 
in general the correlation functions are not continuous at singular configurations in A s [k] 
- indeed, if a sequence of configurations in general positions converges to some singular 
configuration, which implies that some Xi asymptotically reaches lj with j ^ i, the limit of 
the respective correlation functions can be easily seen to coincide with the appropriate one- 
sided correlation function for the limit singular configuration, provided the convergence of 
Xi takes place on one side of lj only, otherwise the limit may fail to exist. The discontinuity 
may therefore arise because the one-sided correlation functions may differ on different sides. 
By Theorem 4 in [18] it follows that for all (l,x) e A 9 [k], k = 1, 2, we have 

<Klx) = l. (15) 

The same paper gives general conditions for this relation to hold for k > 2. Here we are 
interested in the general set-up for k > 2 where it often happens that <f>(l,x) ^ 1. 

Representation theorem for edge correlation functions The following theorem is 
the main result of this paper. 

Theorem 2 For each edge marker collection (I, x) with x C D we have 

<P M {l,x) =Eciop(W^[{l,x)]). 

It is useful at this point to compare Theorem 2 with the simple observation that whenever 
all Xj's in (l,x) belong to a small ball B 2 (xo, r) for some xo G M. 2 then 

<!>(l,x) = N(l,x)(l + 0(r)), (16) 

with N(l,x) defined as in Subsection 4.3, which readily follows by (2,3) and the definitions 
(13,14) of edge correlation functions, taking additionally into account the local boundedness 
of the density m(-) of the activity measure M. with respect to the Haar-Lebesgue measure 
/j, which ensures that the Boltzmann factors exp(— L M (-)) are 1 + 0(r) on B 2 (xo,r), see 
also (12) in [18]. It is interesting to note at this point that, by the consistency of A M , the 
same would hold if we defined N(x, I) to be the cardinality of T(x, I) D D for any convex 
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open D D x and thus N(x, I) does not depend on the field domain D as long as D D x ! 
To proceed, observe that under the same conditions our Theorem 2 yields 

4(1, x) =crop(W° P [(r,x)])(l + 0(r)) (17) 

with W^r[(/, x)\ defined as in Subsection 4.3 and where D r is an r-dependent bounded 
convex domain of diameter 0(r) containing B 2 (xo,r). This is because 0(r) is the proba- 
bility that at least one turning/branching event occurs in the course of the polygonal web 
generating dynamics (W) confined to D r , again in view of the local boundedness of the 
density m(-). To proceed we note that both N(l,x) and crop(W^r [(/, x)}) are, by their 
definitions, invariant with respect to non-singular affine transforms of Consequently, 
upon an appropriate re-scaling we can take r in (16) and (17) arbitrarily small. This way, 
upon comparing (16) and (17) we have established 

Lemma 2 For each edge marker collection (l,x) we have 

cTop(W D [(l,x)]) = N(lx) 
for each D D x. In particular, ciop(yVl,[(l, x)]) does not depend on D D x. 

Note at this point that the relation (17) can be further extended to produce a small r 
expansion, with the coefficient at r k corresponding to instances of polygonal webs with 
exactly k turns/branchings. We do not pursue this topic here though because we are not 
aware of any natural geometric interpretations for the higher order terms of this expansion 
in the style of Lemma 2. 

Relaxing the stop if separated rule As has already been remarked above, in contrast to 
the remaining rather strict dynamic rules, the (W:StopIfSeparated) rule can be some- 
what relaxed without affecting the validity of our Theorem 2. This is made more specific 
in the remark below. 

Remark 1 A direct inspection of the proof of our representation Theorem 2 below shows 
that the result stays valid if the (W:StopIfSeparated) rule, requiring the web branch 
growth to cease immediately when its tip separates from the remaining ones, get replaced 
some other rule where the growth is stopped only on separation but not necessarily imme- 
diately at separation. The only natural constraints are that 
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• At each time moment in the course of the graphical construction the decision on 
whether to stop the growth of a separated branch either depends deterministically 
on the present configuration of branches or at least it is independent of the future 
evolutions of branches given the current branch configuration. 

• Each branch eventually dies before or at the moment s of becoming tangential to the 
current domain boundary dD s . 

The first condition precludes unwanted dependencies whereas the second one is indispensible 
for the technical correctness of our constructions (the above tangency point is the point 
where the time flow direction changes along a branch, and the growth only occurs forward 
in time in our constructions). 

Note that a particular simple example of a stopping rule satisfying the above condi- 
tions is to kill each branch exactly at the time s when it becomes tangential to the current 
boundary dD s . 

In fact, an even deeper analysis of our argument below shows that further relaxation of 
the considered rule are admissible. We do not discuss these details in this paper though 
as they are of no use for our present purposes. 

6 Proof of the representation Theorem 2 

The purpose of this section is to prove our main Theorem 2. Our argument splits into 
several parts and requires some additional concepts. 

Edge marker process To proceed towards establishing our representation theorem for 
edge correlation functions, we shall introduce a Markovian edge marker process whose 
construction can be to some extent regarded as a backwards version of the dynamic repre- 
sentation discussed in Section 3. Roughly speaking, the dynamic representation involved 
an explosion of the field from a single point up to the entire domain D, whereas the edge 
marker process represents an evolution back in time and thus an implosion of the marker 
configuration, eventually to reach one of possible null states. In fact, the edge marker pro- 
cess will be seen to encode the construction of the polygonal web, see (18) below, providing 
an interpolation between the original marker configuration (/, x) and the full polygonal web 
whence the backwards time flow direction. To make all this specific, take the 
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increasing family (D t )te[o,i] of convex compacts satisfying (Dl-4) as chosen in Section 4. 
Next, consider a continuous time branching edge marker process \P S := \1/ s; d, s G [0,1], 
taking its values in finite families of signed and possibly empty edge marker configurations, 
with generic notation 

* s = U p) ■■ (J {p \s),xW(s)) = (t\s),x ( i p \s))t 1 } m 

K. ) p=l 

with m and k p allowed to depend on the time s and where 77^ G {+1,-1}. As the 
notation suggests, the signs 77^ are attributed once and for all to their respective marker 
configurations and do not evolve in time. In addition, we always require that 

(EM:DomainShrink) for each s G [0, 1] the marker points x[ p \s) are all contained in 
the set Di- S , 

that is to say the domain of the process shrinks over time as informally discussed above. 
Given the initial state ^ with all marker points contained in D = D 1 and, in addtion, 
assumed not to contain three different edge marker lines meeting at one point, the process 
^ s is gouverned by the following Markovian dynamics (EM), clearly preserving the latter 
property in view of the absolute continuity M. <C //. 

(EM:DiscardIfSeparated) If at some time s an edge marker (Z-^(s), x\ p \s)) has the 
property that Z-^(s) does not hit the convex hull 

C^ p \s) := conv({x<f\s), q = 1, . . . , m;j = 1, . . . , k q } \ {xf(s)}) 

generated by all the remaining x ( J 9 \s) , s in all marker configurations (l^(s), x^(s)), 
q = 1, . . . ,m, in which case we say that the marker (7^(s), x^\s)) separates from 
\l/ s , then remove the marker (ll P \s) , x\ p \s)) from its configuration (l^ (s) , x^ (s)) . 

If the removal of (l[ p \s) , xf'^s)) makes some other markers separate from \l/ s , the 
(EM:DiscardIfSeparated) rule applies for them as well and they are subsequently 
removed. 

(EM:FoldInwards) Between the critical moments listed below, at each time s > each 
edge marker point x[ p \s) lying at the boundary dD s - in which case we say the 
marker is in boundary position - moves along its corresponding marker line l^ p \s) so 
as to always stay at the boundary of the shrinking domain, that is to say x t f\s + ds) 
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arises as the intersection of if '(s) with dD s+ds . Note that marker points not lying 
at the boundary dD s do not move until they are eventually met by the boundary at 
some later time. 

(EM:DiscardOnCollision) If at some time s two edge marker points x ( f\s) and x^\s) 
within the same marker collection (l^(s), x^(s)) collide (meet) along non-colinear 
directions l!f\s) and ^(s), the two markers are removed from the configuration. 

As a result of the above (EM:DiscardOnCollision) rule some other edge markers 
may separate from \1> S - in such a case the (EM:DiscardIfSeparated) rule applies 
immediately and the markers get discarded. 

(EM:Kill) In the course of the time interval [s, s + ds], an edge marker (l^\s), x ( f\s)) in 
boundary position, moving along the segment x*f\s, s+ds] := xf'\s)x ( f\s + ds), gets 
removed from its configuration (¥ p \s), x&\s)) with probability A^([[x^[s, s + ds]]]). 
These updates are performed independently for all different boundary edge mark- 
ers throughout all marker configurations constituting the process \I/ S , yet they are 
performed simultaneously for all equal edge markers contained in different configura- 
tions, that is to say if (x ( f\s), l!f\s)) = (x^\s), lj*\s)) for some p ^ q then the kill 
events during [s, s + ds] coincide for both these markers. In other words, the killing 
mechanism is a.s. identical for all instances of an edge marker present in different 
configurations constituting \J/ S . 

As a result of the above (EM: Kill) rule some other edge markers may separate from 
^ s - in such a case the (EM:DiscardIfSeparated) rule applies immediately and 
the markers get discarded. 

(EM:TurnAndBranch) In the course of the time interval [s, s + ds], for each boundary 
edge marker (l[ p \s) , x ( f\s)) moving along the corresponding segment x^[s, s + ds] := 
xf'\s)x < f\s + ds), with probability M.(dl) for I e [[x^[s, s + ds]]] a turn-and-branch 
update occurs in the direction of I, which results in replacing the original marker 
configuration (/^,x^) by three offspring marker configurations in which are: 

• : (1^ (s) , x^ (s)) (unmodified offspring), 

• 77W : (i»(s),iW(s)) \ {(l\ P \s),x^\s))} U {(/ t (p) (s + ds) := l,x ( f\s + ds))} (di- 
rectional update offspring - the original marker line l^f\s) turns in the direction 
of I, that is to say ^{s + ds) := I), 
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• -r]^ : (&\s),x (p \s))U{(l,x ( i p \s+ds))} (branched offspring - both the original 
marker and its directional update are present). 

As in the case of (EM: Kill) above, these updates are performed independently for all 
different boundary edge markers throughout all marker configurations constituting 
the process yet they are performed simultaneously for all equal edge markers con- 
tained in different configurations, that is to say if (x ( f\s),l ( f\s)) = (x^\s), lf\s)) 
for some p ^ q then the turning/branching updates during [s, s + ds] coincide for both 
these markers. In other words, the turning/branching mechanism is a.s. identical for 
all instances of an edge marker present in different configurations constituting ^/ s . 

Whenever the above (EM:TurnAndBranch) update is performed with lC\C\ p \s + 
ds) = 0, which results in the directionally updated edge marker separating from 
\l/ s , the rule (EM:DiscardIfSeparated) applies immediately to the corresponding 
directional update and branched offsprings. 

(EM:ForcedTurnAndBranch) If a boundary edge marker point x ( f\s) crosses some 
edge marker line lj P \s), j ^ i, lf\s) ^ ^(s), during the time interval [s,s + 
ds], a forced turn- and-branch update occurs in the direction of ^(s), which results 
in replacing the original marker configuration (1^,2^) by three offspring marker 
configurations in ty s , which are: 

• : (ftP)(s), x^(s)) (unmodified offspring), 

. rj^ : (ZW( a ),x«( S ))\ {(^(s), x[ p) (s))} U {(l^ (s + ds) := lf\s) , xf> (s + ds))} 
(directional update offspring - the original marker line l^\s) turns in the direc- 
tion of Zj p) (s), that is to say l^\s + ds) := lf\s)), 

• : (J(p\s),x(p\s)) U {(lf ) (s),x\ p) (s + ds))} (branched offspring - both the 
original marker and its directional update are present). 

It should be noted at this point that, unlike in the usual (EM:TurnAndBranch) 
discussed above, here (EM:DiscardIfSeparated) has no chance of becoming appli- 
cable directly upon the update because along /^(s) there always exists a direction 
pointing at C^ p \s+ds), namely that towards x ( j'\s+ds). Moreover, observe also that, 
in effect of the so-defined (EM:ForcedTurnAndBranch) update, x ( f\s + ds) and 
Xj P \s + ds) become coupled both in the directional update and branched offsprings. 
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(EM:UnbreakableCouplings) If at some time s in the course of their evolution two 
marker points x ( f\s) and x^\s) are coupled in their configuration then whenever 
they cease to be so in the original configuration or any of its offspring configurations, 
the coupling-breaker configuration is instantly removed from ty s . Note that this is 
equivalent to the rejection of configurations where 

• a coupled edge marker gets killed in a collision (EM:DiscardOnCollision) or 
in (EM:Kill), 

• a coupled edge marker modifies its direction in directional update offsprings 
arising in (EM:TurnAndBranch,ForcedTurnAndBranch). 

The unmodified and branched offsprings do not break couplings and neither can 
a coupling be broken in (EM:DiscardIfSeparated) because coupled markers are 
never separated since they always point at their pair. Note that in contrast to 
(EM:DiscardIfSeparated) rule, where we discard individual edge markers, here 
we remove entire configurations. 

If at some time moment a marker point x^\s) reaches its coupled x ( f\s), both 
markers coalesce and evolve henceforth as one, in particular all coupling restric- 
tions cease to apply. Such meeting and coalescence may occur at the tangency 
point of the respective directional line to the current domain dD s in which case 
the resulting single marker is instantly discarded from the system by application of 
(EM:DiscardIfSeparated). 

The above construction of the edge marker process may seem rather bizarre at the first 
look, but this is in fact a rather simple object. The moving boundary of the shrinking 
domain D s drives inwards polygonal branching random walks of constituent edge markers. 
The directional updating and branching mechanisms of these walks are determined by the 
activity measure M.. The directions of the walks are always chosen to point at the convex 
hull generated by the remaining marker points in the process. If such a choice becomes 
impossible due to edge marker separation, the marker is discarded. Colliding edge markers 
are also discarded. One further rule is unbreakability of once established marker couplings, 
which is ensured by rejecting coupling-breaker configurations. It is important to note that 
with probability 1, in the course of the dynamics all markers eventually 

• either separate and are discarded in (EM:DiscardIfSeparated), 
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• or disappear in collisions (EM:DiscardOnCollision), 

• or are killed in (EM:Kill), 

• or finally they have their configurations annihilated due to coupling breaks, as an 
application of (EM:UnbreakableCouplings). 

Thus, at time 1 the process consists a.s. of signed empty marker configurations. 

Another crucial observation, readily verified by comparing the (EM) and (W) dynam- 
ics, is that, on the event = {+1 : (l,x)}}, the union of trajectories traced by the 
constituent marker points x ( f\s), s G [0, 1], of \I/ coincides with the web W[(l, x)], that is 
to say 

vO/,s)]= U U U k (p) ( s )>- ( 18 ) 

s€[0,l] (rW(s),s(f)(s))6* s ( x \P\ s ),l^\ s ))e(l(p)(s),x( p \s)) 

Moreover, again by the construction, (the history of) each marker configuration present 
in ^/i bijectively corresponds to a complete, minimal and normal collection of branches 
of W[(l, x)], whereas the complete abnormal branch collections correspond to marker con- 
figurations rejected in (EM:UnbreakableCouplings) (this latter correspondence also 
becomes a bijection as soon as the minimality of branch collections is assumed). 

Recall now our assumption made above stating that the killing, directional updating 
and branching mechanisms, while independent for different markers, do coincide for equal 
markers. This assumption is clearly crucial for (18) above to hold, but could be easily 
lifted without affecting the validity of a significant part of the theory presented below. In 
fact, these mechanisms can also be coupled in any other way as soon as the Markovian 
property of the dynamics is preserved. We do not discuss this issue here though as the 
imposed coupling seems to be the most natural one and leading to simplest formulations. 

Some concern may be raised by the branching nature of the (EM) evolution - a natural 
question is whether no cardinality explosions occur for \l/ s . This possibility is easily excluded 
though, as stated below. 

Lemma 3 For each bounded open convex set D and initial condition there exists 
c[D; < +oo such that 

V se[0 ,i] Ecard(^ s ) < c[D;* ]. 
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To see it use first the relation (18) and the discussion following it to conclude that, for all 
s e [0, 1], the expectation Ecard(^ s ) bounded by 2 number of branches of w[(r,s)] > which is 

the maximum number of possible branch collections. Now, the expectation of this number 
is finite because W[(Z, x)] arises from a (critical) binary branching process evolving during 
a finite time interval. Note that the forced branchings do not cause trouble here because 
they only allow to extend already existing lines born at time s to (at the furthest) the 
opposite side of Di_ s , whereas new lines are only born due to the usual critical branching. 

Correlation process Having constructed the edge marker branching process \l/ s , we are 
now going to compose it with the correlation function to obtain the edge correlation process 
$ s , s e [0, 1]. We put 

V ip) (sW ip) (s),x {p) (s)), se[0,l], (19) 

(i"(p)( s ),x(p)(s))e* s 

which means defining the correlation process to be the sum of correlation functions for all 
marker configurations in ty s taken with their corresponding signs. For formal correctness it 
is convenient to adopt at this point the convention that whenever at some time s the marker 
point x^\s) is in a singular position lying on some lj P \s), j ^ i, then the correlation 
function (j>(l^(s), x^(s)) is interpreted as the one-sided correlation function in which the 
marker line Z^(s) is replaced by half-line [^^(s) indicating the direction where x ( f\s) 
came from just before hitting /^(s). In view of (2,3) and since M. <C fi, with probability 1 
this is equivalent to putting in such case <f>(l^(s), x^(s)) := lim n ^ s _ (f>(l^(u), x^(u)) for 
each s where a singularity is reached. 

Assume now that ^ = {+1 : (l,x)} and recall from the discussion following the 
definition of the (EM) dynamics that by the time 1 the marker process \& reaches a ter- 
minal state consisting entirely of signed empty marker configurations. In view of (18) 
and the discussion following it, each such empty marker configuration has its history en- 
coded by some complete normal (minimal) branch collection in W[(/,x)], whereas com- 
plete abnormal (minimal) branch collections correspond to marker configurations rejected 
in (EM:UnbreakableCouplings). Moreover, the sign rj assigned to each empty marker 
configuration in can be readily verified, by induction in the number of branchings in 
the course of crop graph creation, to be (—1) num t , er of branchings _ Observing that for a 
complete branch collection, the number of branchings is simply the difference between the 
number of branches and number of roots (the latter coinciding with the cardinality of the 
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initial marker collection (l,x)) and recalling that 0(0) = 1 we finally conclude from (19) 
and (8) that 

$ = (f>(l,x), $i = crop(>V[(Z,x)]). (20) 

Martingale property of the correlation process With the notation introduced above, 
we claim that the edge correlation process is actually a martingale. 

Lemma 4 The correlation process (& s )se[o,i] ^ s a martingale with respect to the filtration 
T s generated by the marker process 

In view of the relation (20) Lemma 4 means we have just constructed a martingale in- 
terpolating between (l,x) and crop(>V[(Z, x)]). This immediately implies the assertion of 
Theorem 2 upon putting \I>o := {+1 : (l,x)}. Thus, it remains to establish the crucial 
Lemma 4. 

Proof of Lemma 4 In view of the Markovian nature of the edge marker process ty s , to 
prove Lemma 4 it is enough to establish the desired martingale property at s = 0. Moreover, 
for simplicity we present our argument for the initial value of the marker process 
consisting of a single marker configuration (l^ (0) , x^\0)) := (l,x) = (Zj,Xj)^ =1 , whence 
the general argument for higher cardinality of the initial condition is readily obtained by 
straightforward repetition for all configurations in \& . We assume without loss of generality 
that X\ lies at the boundary dD 1 , for only boundary markers undergo evolution under 
(EM) dynamics. The remaining marker points x iy i = 2, . . . ,k may lie both on dD 1 and 
in the interior of D\. 

Using Lemma 3 in [18] and the definition (13) of correlation functions we see that if 
the line l\ does not hit the convex hull conv({a;2, . . . , x k }) then 

0(/i, xi, l 2 , x 2 ; . . . ; l k , x k ) = 0(/ 2 , x 2 ; . . . ; l k , x k ) (21) 

which justifies the (EM:DiscardIfSeparated) rule. Thus, below with no loss of generality 
we constrain ourselves to the case where the (EM:DiscardIfSeparated) rule does not 
apply during the period [0, ds] of the (EM) evolution. 

To proceed with our argument, we shall use the generalised dynamic representation 
described in Section 3, with the same increasing family of convex compacts (D t ) te [ 0;1 ] as 
that used in the construction of the edge marker process. We also recall that D\ — D where 
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D is the field domain. It should be recalled at this point that the generic time parameter 
t of the dynamic representation is related by t = 1 — s to the usual time parameter s of 
the edge marker and correlation processes. As already signalled above, the idea below is 
to relate the dynamics of the edge marker process with the graphical representation under 
inverted (backward) time flow. 

If there is another marker point x iy i ^ 1, in boundary position with the property that 
l\ and U meet in Di \ Di_ ds then, putting % = 2 for notational clarity, by (13) and the 
(GE:VertexBirth) dynamic rule we have, 

<j>(h, xi, Z 2 , x 2 ; k, x 3 ;... ; Z fc , x k ) = (1 + o(l))0(/ 3 , x 3 ; . . . ; l k , x k ). (22) 

Clearly, the above event corresponds to marker point collision under (EM) dynamics and 
thus (22) justifies the (EM:DiscardOnCollision) rule. Keeping this in mind, below we 
only consider the case where X\ does not collide with other marker points during [0,oLs]. 
To proceed, write the correlation-defining event 

E\dl x ,x x \ ■ ■■ ; dl k ,x k ] := {Vf =1 3 eeEdges (^) x { e e, l[e\ e dh} 

as the intersection of £[dli, x±\ and S[dl 2 , x 2 ; ■ ■ ■ ; dl k , x k ] where 

£[dh,xi] := {3 eeEdgcs( _4Ai) xi G e, l[e] G dh} 

and 

£[dl 2 ,x 2 ; ... ; dl k ,x k ] := {VjL 2 3 ee EdgeB(^) x i e e, l[e\ G dU}. 

We shall also denote by xq the intersection point of the marker line l\ with dDi_ds- 
With this notation and taking into account that — A H D arises in the dynamic 
construction with (D t ) te [o,i] as discussed above, we are now in a position to represent 
£[dli,xi, . . . ; dl k , Xk] as the union of the following events, disjoint modulo a set of negligi- 
ble probability, whose names were chosen to represent what happens if we move along the 
field edge covering the marker (li,Xi) under the (EM) dynamics inwards Di and towards 
x . 

(E:GoStraight) £[dl 2 ,x 2 ; . . . ;dl k ,x k ] occurs and an edge ei along l± covering both x 1 
and Xq is present in the field A^. 

(E:TurnOutwards) £[dl 2 ,x 2 ; . . . ;dl k ,x k ] occurs, an edge e\ along l\ covering x\ is 
present but it does not reach xq, instead it turns at some point x' of x\x§ into 
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another field edge e' along a line I' and in the direction consistent with that deter- 
mined by the growth of (D t ), that is to say x' is the first point of e' to be revealed 
by the growing window (D t ). Note that such an edge e' can have only infinitesimally 
small length within D\ — D since when moving from x' in the direction indicated 
by the growth of (D t ) we almost immediately encounter the boundary dD. Thus, we 
can say that e' points outside the domain D and away from Di_ ds , whence the term 
outward turn. Often in such situations we shall also say that e' extends outwards 
from x' along V . 

In terms of the dynamic representation the occurrence of (E:TurnOutwards) is 
equivalent to the occurrence of a (GE:VertexBirth) vertex birth event at x' between 
dynamic representation times 1 — ds and 1, giving rise to the edges t\ along l\ and 
e' along V . 

(E:TurnInwards) S[dl 2 , x 2 ; ■ ■ ■ ; dl k , Xk] occurs, an edge e\ along l\ covering x\ is present 
but it does not reach xq, instead it turns at some point x' of x\Xq into another 
field edge e 1 along a line /' and in the direction opposite to that determined by 
the growth of {D t ), that is to say x' is the last point of e' to be revealed by the 
growing window (D t ). In contrast to the above outward turn, here we turn in the 
opposite inward direction. In terms of the dynamic representation the occurrence of 
(E:TurnInwards) is equivalent to the occurrence of a (GE:DirectionalUpdate) 
at x' where e' extending along /' turns into e± along l±. 

For our considerations below it is convenient at this point to denote by x" the inter- 
section point of the inwards half-line e /_> (starting at x' and determined by e') and 
dDi_ ds . We say that x" lies on I' inwards from x' in such cases. Clearly, x" G e' with 
probability 1 — o(l). 

Observe now that, by the dynamic representation, 

• For (E:GoStraight) we only consider the case that x\ does not collide with any 
other marker point as the converse case has already been handled in (22). In this 
situation we have 

P(E : GoStraight) = (1 + o(l)) [(1 - M([[x^]]))a M [dl 1: x ; dl 2 , x 2 ; . . . ; dl k , x k \ 

(23) 
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- / a M [dl', x"; dh, x ; dl 2 , x 2 ; . . . ; dl k , x k ] 

Ji'e[[xixo]] 

^2 <T M [dl j ,x";dl 1 ,x ;dl2,x 2 ;...;dl k ,x k \\ 

tj-e[[xixo]], Ij^h 

where x" and x" in the above integrals stand for the respective intersection points 
of I' and lj with <9_Di_d s , by definition lying inwards from {x'} :— I' fl x Xi and 
{x'j} := /j fl a^i- To establish (23) note first that on (E:GoStraight) an edge e 1 
in the direction of l\ is present at x at the time 1 — ds of the dynamic construction 
and, in the course of the dynamic construction, during the time period [1 — ds, 1] 
the edge unfolds along x Xi eventually reaching x\. Consequently, there can be no 
(GE:DirectionalUpdate) directional updates along x Xx and neither can there be 
(GE: Collision) collisions with other already existing edges. However, according 
to the (GE) dynamics such directional updates are possible along each line from 
[[xoXi]], consequently the probability that none of these possible turns occurs yields 
the prefactor exp(— «M([[xiXo]])) = (1 + o(l))[l — A^([[xix ]])] whereas the remain- 
ing factor a [dli, Xq; dl 2 ,x 2 ; ■ ■ ■ ; dl k ,x k ] in the first term in the RHS of (23) is the 
probability that the marker points xo,x 2 , ... , x k are covered by their corresponding 
edges as required. We claim that the resulting product 

(1 + o(l))(l - M([[x^]]))a M [dl u x ; dl 2 , x 2 ; . . . ; dl k , x k \ (24) 

represents the probability that the event S[dl 1 ,x ;dl 2 ,x 2 ; . . . ;dl k ,x k ] holds and no 
(GE:DirectionalUpdate) turns occur along x~ox~[. To see it we observe that a pos- 
sible directional update of this kind would yield, during the period [1 — ds, 1] of the 
graphical construction, an outward edge e' of infinitesimal length, almost immedi- 
ately hitting the boundary dDi. Since for almost all time moments in the (EM) 
dynamics the distance between x\ and other marker points x iy i > 1, is strictly pos- 
itive, during its short evolution under the graphical construction dynamics the edge 
e' 

— is overwhelmingly unlikely to have its birth event along the infinitesimal segment 
Xq, x{ affected by the occurrence of E [dl±, x ; dl 2 , x 2 ; . . . ; dl k , x k ], 

— has only a negligible chance of affecting the occurrence of the considered event 
£ [dli, x ; dl 2 , x 2 ; . . . ; dl k , x k ] because, in the course of the graphical construc- 
tion, e' is born after x and other x^s contained in D\-ds get covered by the 
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field and, in addition, e' is, with overwhelming probability, too short to reach 
neighbourhood of any of the remaining x^'s contained in D\ \ D^-ds- 

This nearly independence justifies taking the above product, as required. Next, we 
have to subtract the probability that £[dli, xq; dl 2 , x 2 ; ■ ■ ■ dlk, x k ] occurs and there are 
no directional updates (GE:DirectionalUpdate) along x^xi but a (GE:Collision) 
collision of the edge unfolding from xq along l\ with another already existing edge 
occurs on x x\. There are two possible sources of such collisions 

— XqXi meets an inward edge e' along some V non-colinear with any of the marker 
lines lj . Then the probability of the considered event is 

(1 + o(l))a M [dl',x"; dh,x ; dl 2 ,x 2 ; . . . ; dl k , x k ] 

with the prefactor (l + o(l)) due to the requirements that there be no directional 
updates along x^xi and that the inward edge e' reaches x", which are negligi- 
bly unlikely to fail over the infinitesimal time interval [0, ds]. This expression 
corresponds to the second term in (23) above. 

— XqXi meets an inward edge e' along some marker line lj. In analogy to the case 
above, here the probability of the considered event is 

(1 + o(l))a M [dlj, x"; dh, x ; dl 2 , x 2 ; . . . ;dl k ,x k }. 

This expression corresponds to the third term in (23) above. 

• If Xi is not coupled with any other marker point Xj then for (E:TurnOutwards) we 

have 

P(E: TurnOutwards) = (1 + o(l))M(dl 1 )M([[xTx^\})a M [dl 2 , x 2 ; . . . ;dl k ,x k }. 

(25) 

Indeed, the probability of the vertex birth (GE:VertexBirth) at x' as required for 
(E:TurnOutwards) is M. (dli)Ai(dl'), with the notation as in the above definition 
of the event. Integrating over I' e [[xi^o]] yields Ai(dli)M([\x^x^]]). Moreover, the 
occurrence of such a (GE:VertexBirth) event at some dynamic construction time 
in [1 — ds, 1], as yielding an infinitesimally short outward edge e', is overwhelmingly 
unlikely to affect or be affected by the occurrence of S[dl 2 ,x 2 ; . . . ;dl k ,x k } for pre- 
cisely the same reasons as those justifying (24) above. This nearly independence 
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allows us to express P(E : TurnOutwards) as the product of M(dli)M(dl') and 
a M [dl 2 ,x 2 ; • • • ; dl k ,x k ] with the extra (1 + o(l)) coming also from the requirement 
that there be no (GE:DirectionalUpdate) turns nor (GE:Collision) collisions 
along x'x\. 

Note that the above conclusions are valid regardless of whether some marker line 
lj, j 7^ 1, does cross xqX\ or not. Indeed, if such lj crosses xqXi then the corresponding 
marker point Xj lies in the inward direction from the intersection point for otherwise 
an (EM:DiscardOnCollision) event would occur in the (EM) dynamics, which 
we assumed not to be the case. Thus, a possible outward turn in the direction of 
lj, j ^ 1, occurring along x Xi cannot yield an edge reaching and affecting the status 
of Xj and the occurrence of £[dl 2 , x 2 ; . . . ;dl k ,x k ]. 

• If X\ is coupled with some other marker point Xj, j ^ 1, along li = lj then for 
(E:TurnOutwards) we have P(E : TurnOutwards) = because the probability 
of obtaining in (GE:VertexBirth) at x' an edge t\ exactly colinear with lj is zero 
since M. <C \i. 

• If X\ is not coupled with any other marker point Xj then for (E:TurnInwards) we 

have 

P(E : Turnlnwards) = (1 + o(l))M(dh) / a M [dl', x"; dl 2 , x 2 ; . . . ;dl k ,x k \. 

Ji'emm]] (26) 

Note that if some lj, j ^ 1, crosses x xi then the integral in the RHS of (26) above 
includes the singular term A4(dli)a M [dlj, x"; dl 2 , x 2 ; . . . ; dlj, Xj\ . . . ; dl k , x k ] corre- 
sponding to the situation where V = lj. To establish (26) observe that the probability 
of our edge e\ along l\ arising due to a (GE:DirectionalUpdate) directional update 
at x' on an inward edge e' along I' and of having £[dl 2 ,x 2 ; ■ ■ ■ ; dl k ,x k ] at the same 
time, is At(dli) (directional update probability) times a [dl l ,x"; dl 2 ,x 2 ; . . . ; dl k ,x k \ 
(probability of £[dl 2 , x 2 ; . . . ; dl k , x k ] holding and of having an edge e' along I' ending at 
x' and thus passing through x" modulo negligible measure set) times (l+o(l)) to take 
into account the requirement that there be no further (GE:DirectionalUpdate) 
turns along x'x\ which is satisfied with overwhelming probability. In analogy to our 
previous considerations for (24), also here taking products of the above probabilities, 
modulo (1 + o(l)), is well justified because the directional update at x' has only a 
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negligible chance of affecting the occurrence of £[dl2, x 2 ; ■ ■ ■ dl k , x k ]. Integrating over 
I' yields now (26) as required. 

• If x\ is coupled with some other marker point Xj, j 7^ 1, along l\ = lj then for 
(E:TurnInwards) we have P(E : Turnlnwards) = because the probability of 
obtaining in (GE:DirectionalUpdate) at x' an edge t\ exactly colinear with lj is 
zero since M <^ fi. 

Putting now the above observations and formulae (21,22,23,25,26) together, using the 
definition of the edge correlations (12) and recalling that x\ is in boundary position as 
assumed, we see that, with the notation introduced in the above discussion 

• If the marker (Zi,£i) separates from (l,x) then 

<f>(h, xi, l 2 , x 2 ; ■ ■ ■ ; l k , x k ) = (f)(l 2 , x 2 ; . . . ; l k , x k ). 

• If there is another marker point Xj, % ^ 1, say % — 2, in boundary position and with 
the property that l\ and l 2 meet in D\ \ Di-^ then 

<f>(l u xi, l 2 , x 2 ; l 3 , x 3 ;... ; l k , x k ) = (1 + o(l))0(/ 3 , x 3 ; . . . ; l k , x k ). 

• Otherwise: 

— If x\ is not coupled with any other marker point Xj, j ^ 1, 

a M (dl 1 ,x 1 ; . . . ;dl k ,x k ) = (l + o(l)) [(1 - M([[xJx^]]))a M [dl 1 ,x ; . . . ;dl k ,x k ] + 

M(dl 1 )M([[xTx^\])a M [dl 2 , x 2 ; . . . ;dl k ,x k }+ 

[ [M(dl l )a M [dl',x";dl 2 ,x 2 ;... ;dl k ,x k \- 

Ji'e[[xTxo]] 

a M [dl', x"; dh, x ; dl 2 , x 2 ; . . . ; dl k , x k }} 
^ a M [dl j ,x";dl 1 ,x ;dl 2 ,x 2 ;...;dl k ,x k \\ 

whence, upon recalling the definition (13) of the edge correlation function </>(•), 
4>(li,xi,... ;l k ,x k ) = (l + o(l)) [(1 - M([\xjxii]]))<i>(li,x ; . . . ;l k ,x k ) + 
M([[xTx^]\)(j)(l 2 ,x 2 ; . . . ;l k ,x k )+ 
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/ [<f)(l',x";l 2 ,x 2 ; . . . ;l k ,x k ) - 0(1' ,x";l 1 ,x ;l 2 ,x 2 ; ■ . . ;l k ,x k )}M{dl')+ 

Jl'e[[xTxo]] 

^2 [<j>{lj, x"; l 2 , x 2 ; . . . ; l k , x k ) - x"; h, x ; l 2 , x 2 ; . . . ; l k , x k )]] 

/jG[[xoxT]] 

where the extra positive term in the last sum comes from separate treatment of 
the case V = lj in (26), see the discussion directly following this display. This 
can be further rewritten as 

0(7i,xi; . . . ;l k ,x k ) = (1 + o(l)) [(1 - 2M([[xTx^]\))(f)(l 1 ,xo; . . . ;l k ,x k ) + 

(27) 

M{[[xTx^\])(j)(l 2 , x 2 ; . . . ;k k ,x k ) + 

/ [0(/i, x ;...) + <j)(l', x"; l 2 , x 2 ;...)- <p{l', x"; k,x ; l 2 , x 2 ; . . . )}M(dl')+ 

Ji'e[[xTxo]} 

^2 W^ x "j^ l 2, x 2 ;...] l k , x k ) - (f>(lj,a/j; h, x ; l 2 , x 2 ; . . . ; l k , x k )}] . 

— Likewise, if x\ is coupled with some other marker point Xj, j ^ 1, then 

<f>(h,xi;... ;l k ,x k ) = (1 + o(l)) [(1 - M([[xix^\]))(j)(l u xo; . . . ;l k ,x k ) - 

(28) 

<f>(l',x"; h,x ; l 2 ,x 2 ] ... ; l k ,x k )M(dl')- 

'e[[xixo]] 

^ <f>(h, x'l\ h, x ; l 2 , x 2 \... ; l k , x k ). 

i»e[[xixo]], h^h 

To complete the above discussion we recall that in the complementary case x 1 £ dDi the 
point marker x 1 would not evolve under the (EM) dynamics. 

To proceed we combine (21,22,27,28) and extend these relations for all the other bound- 
ary marker points Xi in (l,x). Recalling the evolution rules (EM) for the edge marker 
process, taking into account that \l/o — {(l,x)} and getting rid of the (1 + o(l)) prefactors 
as ds — > 0, we finally obtain the relation 

$ = (f)(lx) = E J2 v ip) (ds)(f)(l (p \ds),x (p \ds)) = E$ ds . (29) 

(T(p)(ds),x( p Hds))ey ds 

Observe in this context that 

• the second line in equation (27) corresponds to (EM:Kill) rule, 
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• the third line there to (EM:TurnAndBranch) rule, 

• and the fourth line to (EM:ForcedTurnAndBranch). 

Likewise, the absence of certain terms in the coupled version (28) of eq. (27) corresponds to 
annihilation of coupling-breaker configurations in (EM:UnbreakableCouplings). Clear- 
ly, the crucial relation (29) admits straightforward extensions for more general initial con- 
ditions and all time moments between and 1, as discussed at the beginning of our proof. 
Thus, (29) implies in particular that 

<P(l (p) (s),x^(s)), sG[0,l], 

where all r/ p ) signs are converted into pluses, is a positive submartingale. Using Lemma 
3 we see it is a uniformly integrable submartingale. Therefore, noting that < $ s , we 
conclude from (29) that ($ s ) s g[o,i] is a martingale, which completes the proof of Lemma 4. 
□ 

Completing the proof of Theorem 2 With Lemma 4 established, we are now in a 

position to use the martingale representation combined with the relation (20) to complete 
the proof of Theorem 2 as discussed next to the statement of Lemma 4 above. □ 

7 Proof of Lemma 1 

Our proof is based on the so-called defective disagreement loop dynamics developed in 
Subsection 6.4 of [18]. Since Lemma 1 is of a purely technical rather than conceptual 
nature and the quite complicated defective diagreement loop dynamics finds no further 
applications in this paper, we decided not to present its several pages long details here, 
referring the reader to results of Subsection 6.4 and Section 7 in [18] instead. To proceed, 
with (I, x) = (li, , l k , x k ) write 

£[(ZZi,Zi; . . . ; dl k ,x k ] : = {V* =1 3 eeEdges (^) x { G e, l[e] G dk} 

for the correlation-defining events. We claim that for Xj)f =1 G A 9 [k] in general position, 
upon fixing li, x\ \ . . . ; h-i, x k -i, the function 

(lk,x k ) i-> f(£[dl k ,x k \\£[dh,xi] . . . ; dl k -i,x k -i})/ M{dl k ) (30) 
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is well defined, continuous and consequently locally bounded. Clearly, this will imply the 
statement of Lemma 1 upon inductive application for (I, x) G A 9 [k]. The existence and local 
boundedness for (T, x) G A s [k] will then follow as well by noting that the edge correlation 
for a singular configuration is bounded above by the sum of all corresponding one-sided 
correlations and by repeating the existence and local boundedness argument given below 
for the case of one-sided correlations, which is a straightforward repetition omitted here to 
avoid unnecessary technicalities. 

To establish our claim for (30) we use the defective disagreement loop dynamics of 
Subsection 6.4 in [18] with directional updating principle induced by an arbitrary growing 
family (D t )te[o,i] as in (Dl-4) satisfying in addition Do = {x\} so that the anchor point 
A(/i) coincides with x\. In analogy to the proof of Theorem 4 in [18], see also Theorem 10 
there, the conditional law of the polygonal field A on the event £[dli, x\\ . . . ; dlk-i, Xk-i] 
is invariant with respect to the following reversible conditional version of the defective 
disagreement loop dynamics, with s standing for the corresponding time parameter in 
which the dynamics unfolds: 

(Create) With intensity M(dl k )ds, on the event ^£[dl k ,x k ] (that is to say if there is 
no field edge along dl k containing x k ) attempt to emit from x k a disagreement path 
with initial creation phase directed along l k . Should the so generated path result 
in a configuration violating £\dl\,x\\ . . . ; dl k -i, Xk-i], discard the update, otherwise 
accept it. 

(Annihilate) With intensity ds, on the event £[dl k ,x k ] (that is to say if there is a 
field edge along dl k containing x k ) attempt to emit from x k a disagreement path 
with initial annihilation phase directed along l k . Should the so generated path result 
in a configuration violating £[dli,x±] . . . ; dlk-i,x k -i], discard the update, otherwise 
accept it. 

Note that there are no update failures arising due to cycle formation along disagreement 
paths in this dynamics, because the chosen directional updating rule comes from a gen- 
eralised dynamic construction, see Section 6.4 in [18] for details. Denote now by 7r create 
the conditional probability of a succesful Create update attempt during time interval 
(s, s + ds) on the event -^£[dl k ,x k ]. Likewise, write vr annihilate for the respective conditional 
probability for Annihilate update on £[dl k) x k }. Clearly, by detailed balance for £[dl k) x k ], 
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we have 

F(£[dlk,x k ]\£[dh,xi; . . . ;dlk-i,x k -i\) = (1 + o(l))7r crea te Annihilate- (31) 

Note that 7r crcatc < M.{dl k )ds, and that 7r anni M ate > cds for some c uniformly positive 
with respect to small local displacements of x k and l k because there is some positive 
probability that the disagreement path initiated by annihilating the edge at x k does not 
hit any Xi, i < k, and thus does not lead to the violation of £[dl\, x\\ . . . ; dl k -i, x k -i\. 
Consequently, we conclude from (31) that the function in (30) is well defined and locally 
bounded on A 9 [A;]. Its required continuity follows also by (31) in view of the assumed 
continuity of the activity measure density m = dAA/d/i. We have thus established the 
desired properties of the conditional correlation in (30) which completes the proof of the 
lemma. □ 
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